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Eternally inflating universes can contain thermalized regions with different values of the cosmological parame- 
ters. In particular, the spectra of density fluctuations should be different, because of the different realizations of 
quantum fluctuations of the inflaton field. I discuss a general method for calculating probability distributions for 
such variable parameters and analyze the density fluctuation spectrum as a specific application. 



1. Introduction 

In this paper I am going to discuss the struc- 
ture of the universe on super-large scales, so large 
that we are never going to observe them. I shall 
argue, however, that this analysis may help us un- 
derstand some features of the universe within the 
observable range. This is based on the work done 
with my student Vitaly Vanchurin at Tufts and 
with Serge Winitzki at Cambridge University. 

Let me begin with a brief introduction to eter- 
nal inflation. As we know, inflation is a nearly 
exponential expansion of the universe, 



a(t) 



(1) 



which is driven by the potential energy of a scalar 
field <p, called the inflaton. a(t) in Eq.((l|) is the 
scale factor and the expansion rate H is deter- 
mined by the inflaton potential V(ip). Inflation 
ends when (f> starts oscillating about the minimum 
of the potential. Its energy is then dumped into 
relativistic particles and is quickly thermalized. 

A remarkable feature of inflation is that generi- 
cally it never ends completely. At any time, there 
are parts of the universe that are still inflating [ 
|l[ H . The reason is that the evolution of ip is in- 
fluenced by quantum fluctuations. This appies in 
particular to the range of <p near the maximum of 
V((p), where the potential is very flat. As a result, 
thcrmalization does not occur everywhere at the 
same time. We can introduce a decay constant 
r such that r = 1/r is the characteristic time it 
takes ip to get from the maximum to the mini- 
mum of the potential. Then the total inflating 



volume in the universe is proportional to 
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The first factor on the right-hand side describes 
the exponential decay of the inflating volume due 
to thermalization, while the second factor de- 
scribes the exponential expansion of the regions 
which still continue to inflate. For flat potentials 
required for successful inflation, we typically have 
r <C 3H, so that Vinf grows exponentially with 
time. The thermalized volume grows at the rate 
dVtherm/dt = TdV in f /dt, and thus V therm also 
grows exponentially. 

Different thermalized regions in such eternally 
inflating universe may have very different proper- 
ties. Here are some examples. 

The potential V(<p) may have several minima 
corresponding to vacua with different physical 
properties. For example, the values of some con- 
stants of Nature (e.g., the electron mass or the 
cosmological constant) or cosmological parame- 
ters (such as the amplitude of density fluctua- 
tions, the baryon to entropy ratio, etc.) could 
be different in the corresponding thermalized re- 
gions. A more interesting possibility is that the 
"constants" are related to some slowly-varying 
fields and take values in a continuous range. For 
example, the inflaton could be a complex field, 
<p = \ip\exp(ix), with a potential having the 
shape of a "deformed Mexican hat" (that is, with 
some x-dependence) . Then different paths that <p 
can take from the top of the potential to the bot- 
tom will result in different magnitudes of density 
fluctuations Sp/p. The amplitude of the fluctua- 



2 



tions will therefore be different in different parts 
of the universe. Another example is a field x ( un ~ 
related to the inflaton) with a self-interaction po- 
tential U(x)- If U(x) is a very slowly varying 
function of x, then it can act as an effective cos- 
mological constant. Quantum fluctuations will 
randomize x during inflation, and observers in 
different parts of the universe will measure differ- 
ent values of U(x) ■ 

Perhaps the most important example is the 
spectrum of cosmological density fluctuations. 
The density fluctuation Sp/p(l) is determined by 
the quantum fluctuation 5(p(l) of the inflaton field 
ip at the time when the corresponding comoving 
scale I crossed the horizon. Different realizations 
of quantum fluctuations Sip (I) result in different 
density fluctuations spectra in widely separated 
parts of the universe. This uncertainty is present 
in all models of inflation. 

In all these examples, we have parameters x 
which we cannot possibly predict with certainty. 
All we can hope to do is to determine the proba- 
bility distribution V(x)- 

An eternally inflating universe is inhabited by 
a huge number of civilizations that will measure 
different values of x- We can define the probabil- 
ity V(x)dx as being proportional to the number 
of observers who will measure x m the interval 
dx [ ■ Now, observers are where galaxies are, 
and thus V(x)dx is proportional to the number 
of galaxies in regions where x takes values in the 
interval dx- We can then write 

V(x) <x F(xMx), (3) 

where F(x)dx is the fraction of volume in ther- 
malized regions with x m the interval dx, and 
v(x) i s the number of galaxies per unit volume 
(as a function of x)- The calculation of v(x) is a 
standard astrophysical problem, and here I shall 
focus on the volume factor F(x)- 

In this discussion I am trying to avoid the word 
"anthropic", because it makes some people very 
upset, but what I want to emphasize is that the 
approach I have just outlined is as quantitative 
and predictive as it can possibly be. Once V(x) 
is calculated, we can predict, for example, that x 
should have a value in a certain range with 95% 
confidence. 



The first attempts to implement this approach 
encountered an unexpected difficulty. It can be 
traced down to the fact that eternal inflation 
never ends, and the number of galaxies in an 
eternally inflating universe is infinite at t — > oo. 
In order to calculate the volume fraction -F(x), 
one therefore has to compare infinities, which is 
an inherently ambiguous procedure. One can in- 
troduce a time cutoff and include only galaxies 
that formed prior to some time t c , with the limit 
t c — » oo at the end. One finds, however, that the 
resulting probability distributions are extremely 
sensitive to the choice of the time coordinate t 
[ 0, l|. Linde, Linde and Mezhlumian [ |(| at- 
tempted to determine the most probable spec- 
trum of density fluctuations using the proper time 
along the worldlines of comoving observers, which 
they regarded as the most natural choice of the 
time coordinate. They found a probability distri- 
bution suggesting that a typical observer could 
find herself at a deep minimum of the density 
field. On the other hand, if one uses the ex- 
pansion factor along the worldlines as the time 
coordinate, one recovers the standard result [Q]. 
Coordinates in general relativity are arbitrary la- 
bels, and such gauge-dependence of the results is, 
of course, an embarrassment. 

The rest of the paper is organized as follows. 
After reviewing the physics of eternal inflation 
in Section 2, I discuss the spacetime structure of 
an eternally inflating universe in Section 3. This 
will help us understand the origin of the gauge- 
dependence problem. The proposed resolution of 
the problem is discussed in Section 4. As a spe- 
cific application, the spectrum of density fluctua- 
tions measured by a typical observer is analyzed 
in section 5. The conclusions are briefly summa- 
rized in Section 6. 

2. Eternal inflation 

The metric of an inflating universe has a locally 
Robertson- Walker form, 

ds 2 = dt 2 - a 2 (i)dx 2 , (4) 

with the expansion rate given by 

a/a w H(ip) = [8ttV{ip)/3} 1/2 . (5) 
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The potential V(<p) is assumed to be a slowly 
varying function of if. As a result, H is a slowly 
varying function of the coordinates, and we have 
an approximately de Sitter space with a horizon 
distance H~ l . The classical slow- roll evolution 
equation for if is 



ipd « -H'(ip)/4TT. 



(6) 



Quantum fluctuations of ip can be represented 
as a random walk with random steps taken inde- 
pendently in separate horizon-size regions, with 
one step per Hubble time H~ x . The rms magni- 
tude of the steps is 



5(p r 



(H/2n) 



(7) 



We do not have a completely satisfactory deriva- 
tion of this stochastic picture in the general case. 
Its main justification is that it reproduces the re- 
sults of quantum field theory in de Sitter space 
for a free scalar field of mass m <C H (that is, the 
two-point function obtained by averaging a clas- 
sical stochastic field coincides with the quantum 
two point function). For flat inflaton potentials, 
the dynamics of ip should be close to that of a 
free field, so one expects the stochastic picture to 
apply with a good accuracy. 

Let us define the distribution F(<p, t)dif as the 
volume occupied by if in the interval dip at time 
t. It satisfies the Fokker-Planck equation [ pi pi 



d t F + d v J = iH a F 1 



where 



-Ld v (H a+2 F) - ^-H a - X H'F. 

57T^ 4-7T 



(8) 



(9) 



The first term of the flux J describes quantum 
"diffusion" of the field <p, while the second term 
corresponds to the classical "drift" described by 
Eq.(||). The parameter a in Eqs.(||),@ repre- 
sents the freedom of time parametrization, with 
the time variable t related to the proper time r 
according to dt = H l ~ a dr. Hence, a = 1 corre- 
sponds to the proper time, t = r, and a = to 
the scale factor time, t = In a. 

A great deal of research has been done on the 
properties of the Fokker-Planck equation (ph and 



on its solutions. To summarize the conclusions, 
there are some good news and some bad news. 
The good news is that the asymptotic form of 
the solutions of (0) is 



F(<f>,t) -» F{cj>). 



,-yt 



(t -> oo). 



(10) 



The overall factor e 7 * drops out in the normal- 
ized distribution, and thus one gets a stationary 
asymptotic distribution for ip. The bad news is 
that F(tp) has a strong dependence on a, so that 
the results are very sensitive to the choice of the 
time coordinate [ . This is a very disturbing con- 
clusion, and one could have thought that there is 
something wrong with the Fokker-Planck equa- 
tion (||). We shall see, however, that there is a 
good physical reason for the gauge-dependence of 
F(ip), so the equation is not to blame. 

3. Spacetime structure of an eternally in- 
flating universe 

To analyze the spacetime structure of an eter- 
nally inflating universe, one can perform numer- 
ical simulations using the stochastic representa- 
tion of quantum fluctuations [ [l2|, Q]. Here, I 
present the results of the simulations performed 
by V. Vanchurin, S. Winitzki and myself [Q. We 
consider a comoving region which has a horizon 
size I — H^ 1 and a homogeneous inflaton field 
if = at the initial moment t = 0. We evolve tp 
in time increments St according to the rule 

6<p(x) = (p c i(x)St + S<p q (x), (11) 

where <p c i is from Eq.(^) and 5(p q (x) is a random 
Gaussian field of zero mean, 

(Scp q (x)) = 0, (12) 

and with a correlation function 

(Sf q (x)Sf q (x'))=^ (HSt)C(r). (13) 

Here, r is the physical distance between the points 
x and x', C(0) = 1, and C{r) rapidly drops to 
zero at r 3> H^ 1 . The correct asymptotic form 
is C(r) cx r -4 [ in the simulations we set 
C{r) = for r > 2H' 1 . The scale factor a(x,i) 
is evolved according to 



Sa(x) = H(x)a(x)5t. 



(14) 
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The spatial distribution of inflating and ther- 
malized regions in a (2 + l)-dimensional simula- 
tion is shown in Fig.l. We used a double- well 
potential for the inflaton, 

V(ip) = V cos 2 (kv), (15) 

where n = n/2ri, and we only consider the range 
— r\ < <p < r\. Inflating regions are white, and the 
two types of thermalized regions corresponding to 
the minima at (p — ±ry are shown with different 
shades of grey. As time goes on, larger and larger 
fraction of the comoving volume gets thermalized, 
so the inflating regions shrink in comoving coordi- 
nates. The inflating regions can be thought of as 
inflating domain walls separating the two types of 
vacua. Since the domain walls cannot disappear 
for topological reasons, it is clear that inflation 
must be eternal in this model [ [l5| . 

The spacetime structure of the universe in these 
simulations is illustrated in Fig. 2. Now, the ver- 
tical axis is time and the horizontal axis is one of 
the spatial directions. The boundaries between 
inflating and thermalized regions, which play the 
role of the big bang for the corresponding ther- 
malized regions, are infinite spacelike surfaces. In 
the figure, these boundaries become nearly verti- 
cal at late times, so that they appear to be time- 
like. The reason is that the horizontal axis in 
Fig. 2 is the comoving distance, with the expan- 
sion of the universe factored out. The physical 
distance is obtained by multiplying by the expan- 
sion factor a(t), which grows exponentially as we 
go up along the time axis. If we used the physical 
distance in the figure, the thermalization bound- 
aries would "open up" and become nearly hori- 
zontal (but then it would be difficult to fit more 
than one thermalized region in the figure). 

Thermalization is followed by a hot radiation 
era and then by a matter-dominated era during 
which luminous galaxies are formed and civiliza- 
tions flourish. All stars eventualy die, and ther- 
malized regions become dark, cold and probably 
not suitable for life. Hence, observers are to be 
found within a layer of finite (temporal) width 
along the thermalization boundaries in Fig. 2. 

It is now easy to see why the probability dis- 
tributions obtained using a cutoff at t = t c are 
so sensitive to the choice of the time variable t. 




(c) ffl 



Figure 1. A two-dimensional simulation for the 
double- well model at four consecutive moments of 
time [ S3). We evolved a comoving region of ini- 
tial size I = H^ 1 with the initial value of <j> = 0. 
Inflating regions are shown white, while thermal- 
ized regions with <p = +77 and <f> = —r\ are shown 
with different shades of grey. Thermalized regions 
of the same type can join in the course of the sim- 
ulation. For example, regions labeled A and B in 
snapshot (c) have merged into a single region in 
snapshot (d). However, regions of different types 
cannot merge: they are separated by inflating do- 
main walls. 



Any spacelike surface S can be an equal-time sur- 
face t — t c with an appropriate choice of t. De- 
pending on one's choice, the surface S may cross 
many thermalized regions of different types (e.g., 
for t = t), may cross only regions of one type, or 
may cross no thermalized regions at all (say, for 
t = ip with tp in the deterministic slow-roll range) . 
These possibilities are illustrated in Fig. 3 by sur- 
faces iSi,iS2 and £3, respectively. If, for exam- 
ple, one uses the surface S 2 as the cutoff surface, 
one would conclude that all observers will see the 
same vacuum with 100% probability. With a suit- 
able choice of the surface, one can get any result 




Figure 2. Spacetime structure in a one- 
dimensional simulation for the double- well model. 
It can be thought of as a spacetime slice through 
the (2 + l)-dimensional simulation illustrated in 
Fig.l. The shading code is the same as in Fig.l. 



for the relative probability of the two minima. 

We have also performed simulations for a two- 
field model with a potential 

V{ip,x) = K)Cos 2 (^)[l+A(l+cos/?x) sin 4 «<p],(16) 

where Vb, K,/3 and A are constants. We start the 
simulation in a horizon-size region with tp = x — 
at the initial moment t = 0. A spacetime slice 
through one of the simulations is given in Fig. 4, 
where different values of x are shown with differ- 
ent shades of grey in the inflating regions, while 
the thermalized regions are left white. The field 
X is homogeneous at the bottom of the figure 
(t = 0), but gets randomized by quantum fluctu- 
ations as time goes on. Different regions of space 
thermalizc with different values of x, and x takes 
all of its values on each thermalization surface. 
It is clear that the simple explanation of gauge- 
dependence that I gave for the one-field model, 
that constants surfaces can be chosen so that 
they cross one type of thermalized regions and 
avoid the other, does not apply in this case. For 
the two-field model, the situation is more subtle. 
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Figure 3. Si is a surface of a constant proper 
time. It crosses many thermalized regions of 
different types. £>2 is a spacelike surface which 
crosses regions of only one type. S3 is a space- 
like surface which does not cross any thermalized 
regions. 



Note that the thermalized regions in Fig. 4 have 
a fractal-like pattern, with larger numbers of 
increasingly narrow "spikes" appearing at later 
times. These spikes represent newly-formed ther- 
malized regions, and although they are rather 
small, their number grows exponentially with 
time, and they actually dominate the thermalized 
volume. A change of the time variable t results 
in a deformation of the cutoff surface, accompa- 
nied by a (substantial) change in the population 
of the newly-formed regions that are being in- 
cluded. This is the origin of the gauge dependence 
of the cutoff procedure. 

4. The proposal 

The resolution of the gauge dependence prob- 
lem that I proposed in Ref. [ [l6| is to calculate 
the probability distribution for x within a single, 
connected thermalized domain. Each thermalized 
domain can be infinitely extended and contains an 
infinite number of galaxies, but it is sufficient to 
use a large finite part of the domain. If the field x 



2n/3 




Figure 4. Spacetime structure in a one- 
dimensional simulation for the two-field model 
showing the evolution of a comoving region of 
initial size I = H^ 1 with initially homogeneous 
= and x = "7 ' P- The values of the field x 
are shown throughout the inflating region using 
different shades of grey. The shading code is in- 
dicated in the bar on the right of the figure. The 
thcrmalized regions are left white. 



varies in a finite range, it will run through all of its 
values many times in a sufficiently large volume. 
We expect, therefore, that the distribution F(\) 
will converge rapidly as the volume is increased. 
It does not matter which thermalized domain we 
choose to calculate probabilities: all domains are 
statistically equivalent, due to the stochastic na- 
ture of quantum fluctuations in eternal inflation. 
This is a very simple prescription, and I am a 
bit embarrassed that I did not think of it earlier, 
having thought about this problem for a number 
of years. 

With this prescription, the volume distribution 
F(x) can be calculated directly from numerical 
simulations, and we have done that in [ |l3| for 
the two-field model (16). In some cases an an- 
alytic calculation is also possible. Suppose, for 
example, that the potential V(ip, \) is essentially 
independent of x f° r W\ < ^Oj where ipo is in 



the deterministic slow-roll range, where quantum 
fluctuations of <f> and \ can be neglected com- 
pared to the classical drift. (The corresponding 
conditions on the parameters of the potential ( |l6| ) 
are specified in [ [l(| Then, the evolution of 
ip at ip > ifo is monotonic, and a natural choice of 
the time variable in this range is t = ip. The prob- 
ability distribution for x on the constant-" time" 
surface ip = (po is 



Fo(x) = F(fo,X) = const, 



(17) 



since all values of \ are equally probable at 
ip < fQ. We are interested in the probabil- 
ity distribution on the thermalization surface, 
F(x) = F(ip*,x), where <p* is the value of ip at 
thermalization. This is given by [ [ll| 



F(x)<xF ( X o)eM3N(xo)]det 



&X 



(18) 



Here, xo is the value of x at V = Vo that classi- 
cally evolves into x a * V*, N(xo) is the number 
of e-foldings along this classical path, exp(37V) 
is the corresponding enhancement of the volume, 
and the last factor is the Jacobian transforming 
from xo to %. In many interesting cases, x does 
not change much during the slow roll. Then, 



F(x) oc exp[3AT( X )]. 



(19) 



In a more general case, when the diffusion of 
X is not negligible at ip > ipo, the distribution 
F(x) can be found by solving the Fokker-Planck 
equation with t = (p in the range tpo < ip < (p* 
and with the initial condition (p"7|). The corre- 
sponding form of the equation was derived in [ 



OF 

dtp 



-d 2 

x 



H 3 F 
2ttH' 



-0 y 



E2 



-F 



\2tiH 
H' 



F.{20) 



We have solved this equation for the two-field 
model ([l6]) with the same parameters that we 
used in numerical simulations and compared the 
resulting probability distribution -F'(x) with the 
distribution obtained directly from the simula- 
tions. We found very good agreement between 
the two (see Ref.[ 0| for details). 
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5. Density fluctuations 

As a specific application of the proposed ap- 
proach, let us consider the spectrum of density 
perturbations in the standard model of inflation 
with a single field tp. The perturbations are de- 
termined by quantum fluctuations Sip; they are 
introduced on each comoving scale at the time 
when that scale crosses the horizon and have a 
gauge-invariant amplitude [ m 



Sp/p = 8wH8<p/H', 



(21) 



where H' = dH/dip. With an rms fluctuation 
(Stp) rms = H/2tt, this gives 



(Sp/p) r 



4H 2 /\H'\ 



(22) 



Fluctuations of ip on different length scales are 
statistically independent and can be treated sep- 
arately. We can therefore concentrate on a single 
scale corresponding to some value tp — tpo, disre- 
garding all of the rest. 

On the equal-" time" surface <p = tpo, the fluc- 
tuations S(p can be regarded as random Gaussian 
variables with a distribution 



Fo(Sip) oc exp 



H. 



o 



(23) 



where Hq — H(<po). We are interested in the 
distribution F(8ip) on the thcrmalization surface 
(p = ip% . This will be different from F if there is 
some correlation between Sip and the amount of 
inflationary expansion in the period between tp 
and ip*. In fact, there is such a correlation. If tp 
fluctuates in the direction opposite to the classical 
roll, then inflation is prolonged and the expansion 
factor is increased. Otherwise, it is decreased, 
and we can write 



F(5ip) oc F (5ip) exp(3£T <ft), 
where 

St = -(4ir/H' Q )5tp 



(24) 



(25) 



is the time delay of the slow roll due to the fluc- 
tuation Sip. 

Combining Eqs.(p3|)-(p5|), we obtain [ p6) 



F(S(p) oc exp 



2^ 
H 2 



(Sip - 5'ipf 



(26) 



which describes Gaussian fluctuations with a 
nonzero mean value, 

6ip = 3H$/nH' . (27) 

This is different from the standard approach [ [t]] 
which disregards the volume enhancement factor 
and uses the distribution (^3|). The effect, how- 
ever, is hopelessly small. Indeed, 



Sip 



(Sip) 7 



6H§ 
H' 



6p> 

P / rms 



1CT 



(28) 



We thus see that the standard results remain es- 
sentially unchanged. 

6. Conclusions 

Eternally inflating universes can contain ther- 
malized regions with different values of the cos- 
mological parameters, which we have denoted 
generically by \. We cannot then predict x with 
certainty and can only find the probability distri- 
bution V(x)- Until recently, it was thought that 
calculation of V inevitably involves comparing in- 
finite volumes, and therefore leads to ambiguities. 
My proposal is to calculate V in a single thermal- 
ized domain. The choice of the domain is unim- 
portant, since all thermalized domains are sta- 
tistically equivalent. This apprach gives unam- 
biguous results. When applied to the spectrum 
of density fluctuations, it recovers the standard 
results with a small correction O(10~ 5 ). 

It should be noted that this approach cannot be 
applied to models where x is a discrete variable 
which takes different values in different thermal- 
ized regions, but is homogeneous within each re- 
gion. [An example of such a model is the one- field 
model (|l5|) .] One can take this as indicating that 
no probability distribution for a discrete variable 
can be meaningfully defined in an eternally in- 
flating universe. Alternatively, one could try to 
introduce some other cutoff prescription to be ap- 
plied specifically in the case of a discrete variable. 
Some possibilities have been discussed in [[It], |l8| |. 
This issue requires further investigation. 
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